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SUMHARY 

IMS  report  discusses  the  problem  of  optima,  determination  of  position 
by  a  navigating  devioe  employing  independent  inertial  and  velocity  meaenre- 
«enta.  The  measurements  are  assumed  to  be  subject  to  random  errors. 

The  problem  is  related  to  the  theory  of  statistical  estimation  as  veil 
as  to  filtering  theory.  Explicit  „«»  methods  of  position  computation  ar, 
derived  for  several  special  cssee.  Formulas  are  derived  by  means  of  which 
the  variance  of  the  error  in  computed  position  can  be  determined.  ihe 

asymptotic  variance,  for  large  time  of  flight,  of  the  error  in  computed 
position  is  discussed. 
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LIST  OF  SYMBOLS 


Note:  A  detailed  explanation  of  certain  parts  of  the  notation  is  given 
on  p„  9. 


Symbol 

x(t) 

v(t) 

b(t) 

x 

0 

X' 

0 

B(t) 

V(t) 

Xo 

X* 

o 

T 

W 

aliPli02>P2» 

Yl,Y2,Y3,Yii 


C 


*1 

0y(5jt) 

A 
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Definition 
vehicxe  position  at  time  t 
vehicle  velocity  at  time  t 
xM(t)  +  x(t) 
initial  position  of  vehicle 
initial  velocity  of  vehicle 
accelerometer  dial  reading 
velocity  dial  reading 
independent  estimate  of  initial  position 
independent  estimate  of  initial  velocity 
elapsed  time  since  beginning  of  flight 
defined  by  Eq.  (IV.3) 

parameters  describing  instrument  error  statistics 

error 

defined  by  Eq,  (m„6) 
defined  by  Eq.  (IV. 6) 

accelerometer  error  autocorrelation  function 
velocity  dial  error  autocorrelation  function 
frequency  of  Earth's  radius  pendulum 
defined  by  Eq.  (III. 3) 
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A.  BACKGROUND 

Considerable  study  has  been  devoted  in  recent  years  to  the  problem  of 
navigating  airborne  vehicles  by  utilizing  acceleration-sensitive  devices 
housed  in  the  vehicle.  In  principle,  if  the  position  and  velocity  with 
respect  to  inertial  space  at  some  initial  instant  t  n  0  were  known  exactly, 
and  if  the  measuring  and  computing  devices  in  the  vehicle  functioned  with¬ 
out  error,  it  would  be  possible  for  a  computer  in  the  vehicle  to  calculate 
the  vehicle's  exact  position  at  any  instant  t  D  T>0.  To  decrease  the 
effect  of  accelerometer  errors,  systems  incorporating  a  combination  of  both 
accelerometers  and  independent  velocity  measuring  devices  (such  as  doppler 
radars)  have  been  proposed,  and  made  the  subject  of  a  well  developed  liter¬ 
ature.*  In  this  report  we  will  focus  our  attention  on  the  accuracy  achiev¬ 
able  with  such  a  combination  type  of  system. 

The  computed  position  contains  an  error  due  to  the  errors  in  the  various 
components  of  the  system — for  example,  errors  in  the  accelerometer,  in  the 
independent  measurement  of  velocity,  in  gyros  or  star  trackers  used  to  main¬ 
tain  a  reference  direction,  in  the  altimeter,  in  initial  values  of  position 
and  velocity,  etc.  The  question  naturally  arises:  In  a  position  computer 
utilizing  independent  accelerometer  and  velocity  measurements  what  is  the 
"best"  way  to  combine  the  information  from  these  sources  to  yield  a  computed 
position?  In  order  to  make  a  solution  of  this  problem  feasible  it  is  necessary 
to  specialize  the  problem  a  good  deal.  Only  a  few  sources  of  error  are 


See  the  bibliography  of  Ref.  2. 
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considered)  various  other  simplifying  assumptions  are  made.  The  problem 
thus  specialized  is  still  of  considerable  interest;  many  of  the  papers 
citea  make  the  same,  similar,  or  more  restrictive  assumptions. 

This  report  will  be  restricted  to  the  presentation  of  the  theoretical 
work.  Numerical  computations  based  on  the  theory,  and  using  experimental 
noise  statistics,  are  being  planned,  and  the  resulting  data  will  be  presented 
in  a  later  report.  One  of  the  objectives  of  the  numerical  work  will  be  to 
compare  the  optimum  system,  which  turns  out  to  have  time-varying  properties, 
with  mure  conventional  non-time-varying  systems, 

D ,  GENERAL  NATURE  OF  INPUT  DATA  TO  THE  COiOTTER 

For  expository  purposes  we  will  first  consider  an  ideal,  error-free 
system,  and  then  approximate  to  the  actual  system  by  adding  random  noise. 

In  the  ideal  system  the  vehicle  in  question  is  assumed  to  move  in  a 
great  circle,  at  constant  altitude  above  a  spherical  and  non-rotating  earth. 
The  task  of  the  computer  is  to  calculate  the  position  x(t)  along  the  great 
circle. 

In  what  follows  capital  letters  will  denote  variables  which  include 
random  error  components,  and  the  corresponding  lower  case  letters  will 
denote  the  same  variables  without  the  error  components.  The  error  compon¬ 
ents  will  be  denoted  by  G  vdth  an  appropriate  subscript. 

We  will  refer  to  the  variables  upon  which  the  course  computer  operates 
as  dial  readings".  Altogether,  there  are  four  dial  readings  (Fig.  1)~ 
acceleror  ter,  velocity,  initial  position,  initial  velocity.  The  first 
two  vary  during  i light;  the  latter  two  are  preset  before  commencement  of 
flight  and  merely  serve  as  a  permanent  store  of  the  preset  information. 
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L 


Accelerometer  dial 
3(t)  -  b(t)  +  £R(t) 


Velocity  dial 


Initial  position 

estimate  , 

X  01  x  +  0., 
o  o  X 


Initial  velocity 

estimate  c 

X  1  ■  x '  +  Cv  , 
o  o  X ' 
o 


Fig.  1 


(a)  Accelerometer  Dial  Reading 

As  used  in  this  report,  the  term  "accelerometer11  is  actually  somewhat 

of  a  misnomer,  as  the  apparatus  in  question  not  only  serves  to  indicate  the 

resultant  force,  but  also  incorporates  a  provision  for  maintaining,  by  means 

of  gyros  and/or  star  trackers,  a  fixed  direction  in  inertial  space.  Theory 
# 

shows  that,  assuming  no  instrumental  errors,  the  accelerometer  dial  can  be 
calibrated  so  that  at  time  t  it  will  read 

(1.1)  b(t)  “  x"(t)  +  _Q2x(t)  . 

In  line  with  the  remar'k  at  the  beginning  of  the  paragraph  the  actual 
reading  B ( t )  will  be  assumed  to  differ  from  b(t)  by  an  additive  noise  term 
63(t): 

d.2)  B(t)  -  b(t)  +L(t)  =  x"(t)  +:  ,2x(t)  +L(t) 

tf  0 

Cf.  e.g.  Ref.  1. 
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Jial  Heading 


A  true  velocity  measuring  device  would  indicate 


v(t)  “  x'(t) 


The  actual  reading,  however,  will  be  of  the  form 


V(t)  *•  v( t)  +  0^( t)  “  x'(t)  +  Oy(t) 


(c)  Initial  Position  Estimate 


Tha  analysis  admits  the  possibility  of  an  estimate  X  of  the  initial 


position  x  being  set  into  the  computer  before  commencement  of  flight. 


("Initial"  refers  to  the  instant  when  the  course  computer  begins  to  function.) 


We  put 


X  B  X  +  Gy 
o  o  A 


where  C  v  is  the  error  of  X 
X  o 

o 


(d)  Initial  Velocity  Estimate 


In  addition  to  X  an  estimate  X1  of  the  initial  velocity  x'  may  be 
0  0  o 


preset  into  the  computer.  As  above,  we  put 


X  1  B  x 1 
o  o  X 1 


4  short  remark  apropos  X  ana  X'  may  serve  for  proper  orientation  of 
r  0  0 


> 
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the  reader  at  this  point.  Every  bit  of  information  available  to  the 

computer  may  be  expected  to  aid  it  in  its  estimate  of  position.  The 

computer  could  estimate  xq  and  just  on  the  basis  of  the  accelerometer 

and  velocity  dial  readings  B(t)  and  V(t)  by  using  the  estimate  V(0)  for 

x  *  and  |B(0)  -  V'(0)  }tct  for  x  ,  X  and  X  1  are  additional  estimates 
o  L  |  /  o  o  o 

which  are  independent  of  the  accelerometer  and  velocity  dial  readings. 

Once  the  accelerometer  and  velocity  dial  readings  become  available,  the 
computer  has  essentially  two  independent  estimates  of  xq  and  of  xj  at  its 
disposal.  The  two  estimates  of  x^  ,  or  of  x^  ,  can  be  combined  (e.g,  by  a 
linear  weighting)  to  obtain  a  new  estimate  which  is  better  than  either  of 
the  original  ones  taken  by  itself.  It  is  thus  clear  that  (except  in  degen¬ 
erate  cases)  the  separate  specification  of  Xq  and  X^  is  not  superfluous, 
but  rather  these  quantities  may  be  expected  to  play  a  role  in  the  optimum 
position  estimation  process. 

C.  ERROR  STATISTICS 

The  accelerometer  and  velocity  dial  errors  are  assumed  to  be  random 
processes,  with  autocorrelation  functions  y^(s,t)  and  0y(s,t)  respectively; 
and  with  zero  mean  over  a  statistical  ensemble  of  flights.  That  is,  for 
every  instant  t, 


(1.7) 


5B(t)  =  0 


£v(t)  «■ o 


and,  for  each  pair  of  instants  (s,t), 

£g(t)£g(s)  "  0g(s, t) 

cy(t )£,y(s)  -  0v(s,t) 
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The  zero  mean  holds  l'or  a  statistical  ensemble  of  flights.  The 

possibility  exists  that  on  any  single  flight,  the  dJ  1  "actings  may  be 

biased.  The  average  value  of  such  a  bias  will,  <''■  t  a  large  number  of 

flights,  be  zero.  The  existence  of  such  biases  causes  a  positive  constant 

term  to  be  added  to  the  autocorrelation  functions. 

The  assumption  of  zero  ensemble  mean  for  the  errors  is  clearly  no 

restriction  on  generality,  since  if  the  mean  were  a  known  quantity  different 

from  zero,  it  could  be  subtracted  from  the  dial  readings. 

The  quantities  L  and  tj,  of  (1.5)  and  (1.6)  are  assumed  to  be  (over 
o  o 

an  ensemble  of  flights )  random  variables  with 

(1.9)  "  0  j  c2  .  r,  J  ^r,  "  0  5 

o  oo.  0 

The  dial  errors  are  all  assumed  to  be  statistically  independent  of  each 
other  and  of  the  true  path. 


D.  THE  OPTIMIZATION  CRITERION 

The  optimization  criterion  is  as  follows:  at  each  instant  T  ,  the 
computer  should  form  an  estimate  of  position  x(T),  based  on  the  dial  readings 
for  0  ■  t  *  T,  such  that 

a)  x(T)  -  :c(T),  the  error  in  position  estimate,  is  independent  of  the 

path.  This  is  equivalent  to  assuming  that  no  a-priori  information  e.g. 
statistical  information — about  the  set  of  possible  paths  is  utilized. 


b)  x(T)  depends  in  a  linear  manner  on  the  dial  readings. 

c)  Expected  value  (with  respect  to  a  statistical  ensemble  of  flights)  of 

•  '  ]2 

'x(T)  -  x(T)  is  a  irdninrum  for  all  estimates  satisfying  (a)  and  (b). 
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S°  fUrtiler  r“trlotl»-"-s  «Jo  on  the  possible  ways  of  coining  ? 
the  information  contained  l„  ths  dial  readings. 

^saaag-gitsisg^omMi^  s„lutlon  t0  th>  g2ag 

(linear  or  nonUng^  -m..-.» 

estimates  of  position. 

Conditions  (,)  and  (b)  together  aith  the  assumed  dial  error  statistics 
inp]y  thai  the  optimum  position  estimate  will  satisfy 

^(T)»  wiieia  °-x(t)  haa  Zero  ensemble  mean,  depends  in  a 
linear  mnrer  on  the  ^  errors,  and  does  not  depend  on  the  path.  (Consequently, 
the  system  is  dynamically  exact.) 

S.  OUTLINE  OF  Trig  REMAINDER  OF  THE  REPORT 

Section  I!  contains  an  exposition  of  the  general  method  of  attack  on 
the  problem  at  hand. 

It  happens  that  detailed  resell,  are  easy  to  obtain  in  tuo  special 
cases,  viz,; 

SffiS-JL  ■  S^g(s,t)  *>  constant  ■>  y^ 


Case  2 


0y(s,t)  “  constant  »  y 


■IWe  can  be  interpreted  as  folio™,  either  the  accelerometer  dial 
(Caae  1)  or  the  velocity  dial  (Case  1)  is  noiseless  except  for  bias. 

section  in  la  devoted  to  the  treatment  of  Case  1,  Section  IV  to 
aae  2.  Explicit  solutions  for  the  optimum  position  estimate  are  obtained 
only  after  oven  further  specialization,  to  wit: 

Explicit  solution  obtained  for 


(1. 10) 


0v(s,t)  »  x2  *  p  e  °2 


3-t 


Over  an  ans  entile  of  flights. 
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d*!!)  0B(s,t)  “  +  Pje  1  '  ' 

Ihe  application  of  our  methods  to  obtain  explicit  optimum  position 
estimates  for  other  autocorrelation  functions  is  further  discussed  in  the 
sections  dealing  with  Cases  1  and  2. 

Sections  III  and  IV  also  contain  expressions  from  which  can  be  obtained 
the  variance  of  the  error  in  'x(T)  as  a  function  of  T. 

The  results  for  Cases  1  and  2  are  summarized  at  the  end  of  Sections  III 
and  IV,  respectively, 

A  third  special  case  which  is  treated  in  Appendix  II  is 

Case  3  .  , 

f  .  .  ~cu  s-t 

S^(3,t)  «  pie  1  ! 

)  -a,  s-t 

/  07(s, t)  -  pge  2 

The  solution  to  this  case  is  relegated  to  an  appendix  because,  although  it 
is  qualitatively  interesting,  it  is  so  complicated  as  to  be  almost  useless 
for  purposes  of  computation.  The  solution  and  the  method  used  to  obtain  it 
are  outlined,  but  not  given  in  full  detail. 

Section  V  shops  that  the  proble.  ot  opt™  poeitlon  eatlmaon  can  be 
reduced  to  a  linear  least  squares  filtering  problem. 

Section  VI  contains  a  discussion  of  asymptotic  errors. 
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(a) 

(b) 


T  is  the  elapsed  time  since  the  beginning  of  the  flight. 


Dial  readings  are  denoted  by  B,  V,  XQ,  X^ 


The  values  these 

quantities  would  have,  if  the  dials  were  errorless,  are  denoted  by 


(c) 


(d) 


(e) 


(f) 


(g) 


b,  v.  x  ,  x1  . 

For  0  »  t  *■  T,  the  computer 'o  estimate  of  x(t)  based  on  all  the 
dial  readings  up  to  time  T  is  denoted  by  x(tjT).  Thus,  x(0;T)  repre¬ 
sents  the  computer's  estimate  of  initial  position  based  on  all  dial 
readings  up  to  time  T.  The  quantity  'x(TjT)— i.e.  the  computer’s 
estimate  of  position  at  time  T  based  on  the  dial  readings  up  to 
time  T— will  be  denoted  byx(T), 

Similar  notation  is  used  for  the  computer  's  estimates  of  other 
quantities  such  as  velocity:  x'(tjT),  x^(OjT),  etc,* 

Errors  in  all  quantities  are  denoted  by  6  with  an  appropriate 

subscript,  as:  Cfi  «  B  -  bj  »  x(T)  -  x(T)j  etc. 

<* 

However,  errors  such  as 
denoted  by  6g(t)  ,  ^A(tjT)  ,  etc. 


,  ^x(tjT)  *  e^c*>  will  usually  be 


For  functions  of  time  F(t),  quantities  such  as  F(t^)  will  some¬ 
times  be  denoted  by  F^, 

The  operator  — *  +  1  Z  is  denoted  by  ^ :  X  x(t)  ■  x"(t)  +  -£”2  x(t)  . 
d  dt 

The  operator  ^  will  sometimes  be  denoted  by  ra  x'(t). 

For  a  list  of  the  other  symbols  used,  see  p,  v. 


In  the  following  pages,  reference  will  often  be  made  to  the  "best" 
estimates  of  x(t),  x'(t),  etc;.,  based  on  all  the  dial  readings  up  to  time 
T,  for  any  t  in  the  range  0  *  t  **  T.  In  such  cases,  "best"  means  best  in 
the  sense  of  the  optimisation  criterion  of  part  D  (p.  6),  with,  for 
position  estimation,  x( ^replaced  by  x(tjT)  and  x(T)  replaced  by  x(t); 
for  velocity  estimation,  x(T)  replaced  by  x»(tjT)  and  x(T)  replaced  by’ 
x'(t);  etc. 
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II,  METHOD  OF  SOLUTION 

One  approach  to  the  problem  of  finding  estimates  having  minimum  error 
variance  is  the  maximum  likelihood  method.  Hie  role  of  this  method  in 
conventional  estimation  problems  is  treated  in  detail  in  many  standard  works. 

The  method  may  be  briefly  suramed  up  as  follows:  Consider  a  vector  with 
random  components  t?  «  (up  ...  ,  u^)  ,  and  suppose  that  the  joint  density 
function  of  in,  ...  ,  depends  on  some  parameters  dp  ...  ,  a^; 0  Let 
a  "  (op  ...  ,  a^).  Let  the  density  function  be  denoted  by  p^U  |  a^)  .  Then, 
provided  the  density  function  satisfies  certain  regularity  conditions,  the 
unbiased  estimates  of  minimum  variance  of  the  parameters  cl,  ...  ,  o^,  based 
on  observed  values  Up  ...  ,  U^,  are  found  by  maximizing  p(Up  ...  ,  Un  Op  ...  : 
with  respect  to  Op  ...  ,  aft.  The  maximizing  values  of  Op  ...  ,  are  the 
required  estimates.  • 

In  the  next  two  sections,  the  maximum  likelihood  method  will  be  applied 
assuming  Gaussian  statistics.  This  will  give  the  estimate  of  minimum  variance 
among  all  unbiased  estimates.  As  was  pointed  out  in  the  introduction,  however, 
the  resulting  estimate  is  also  the  linear  estimate  of  minimum  variance,  no 
matter  what  the  noise  statistics  are,  provided  the  means  and  correlation 
functions  remain  the  same. 

Another  way  of  expressing  this  is  as  follows:  for  any  statistics,  the 
linear  estimate  of  minimum  variance  can  be  found  by  replacing  the  actual 
statistics  by  Gaussian  statistics  having  the  same  means  and  correlation 
functions,  and  then  applying  the  maximum  likelinood  method.  If  the  statis¬ 
tics  actually  are  Gaussian,  this  gives  the  minimum  variance  estimate  among 
all  unbiased  estimates. 


Cf .  e .g.  Ref.  3. 
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The  application  of  the  maximum  likelihood  method  to  the  problem  at 
hand  is  simpler  to  explain  first  for  Cases  1  and  2, 

In  Case  1,  for  example,  (^(s,t)  -  Yp  so  that  the  accelerometer  dial 
is  perfect  except  for  bias.  This  means  that  the  accelerometer  error  is 
constant  for  each  flight  and  equal  to  the  initial  accelerometer  error. 

Now  suppose  that  the  velocity  dial  readings  are  known  only  at  times 

0  "  to<tl<,,,<tn-l  01  T’  where  ti  ^  equally  spaced  over  the  inteival 

(0,T),  The  vector  1?  “  (B  .  V  V  \r  y  , 

o’  o’  1’  **•  ’  vn_i>  x0>  has  a  density 

function  which,  for  Gaussian  statistics,  can  be  explicitly  written  down  in 
terms  of  0y(s,t),  and  can  be  shown  (Section  III)  to  depend  on  just  three 
unknown  parameters,  which  can  be  taken  to  be  x  ,  x-  and  Xx  -  x"  *-02  v 

O'  O'  O  O  ~  0* 

Also  x(T)  can  be  expressed  in  t.r„  of  xo,  x.,  end  know  factions  of 

T. 

Therefore,  one  must  maximize  p(tf|x  x.,  &  )  vrtth  respect  t 
y ■'  v  v  o'  o’ 

and^xo  (with  the  observed  dial  readings  as  components  of  o). 


The  resulting  estimate 


x(T) 


would  be  the  unbiased  estimate  with 


mini—  error  variance,  if  the  velocity  dial  readings  were  known  only  for 
ti»,s  t..  The  lilt,  l£jj(r;]n  is  the  required  estivate  of  x(T)  hated 
on  all  information  contained  in  the  dial  readings  for  0  “  t  ■  T. 

In  Case  2,  the  velocity  dial  is  perfect  except  for  bias.  In  this  case, 
the  random  vector  under  consideration  is  "u  *  (Bo .  B] ,  . . .  ,  B  V  X  X1 
The  density  function  depends  on  two  unknown  parameters,  xq  and  x^;  x(T)  can 
be  expressed  in  ter™  of  those  and  known  functions  of  T.  IhcrcfL,  one  m. 
maxindEe  p(U  |xq,  xp  with  respect  to  xq  end  so  obtai 
mi  t  ss  n  oo # 


5  and  take  th 


In  the  general  case,  the  problem  can  in  principle  be  approached  as 
follows,  The  four  dial  readings  B(t),  V(t),  X<),  x<;  are  shown  in  section  V 


CONFIDENTIAL 


in  the  sense  that,  knowing  R(t)  and  Q(t)  over  some  interval  0  “  t  T, 

one  can  derive  V(t),  B(t),  Xq,  X^  for  that  interval;  and  vice  versa. 

The  readings  of  the  Q  and  R  dials  at  times  0  •»  t  ^  t,  <t  ,  ®  T 

ox  n-1 

form  a  2n-dimensional  random  vector,  the  density  function  of  which  depends 

on  the  vector  ~x  c  (xqI  ...  ,  xn_^);  here  *»  x(t^).  The  autocorrelations 

and  cross  correlation  of  Q(t)  and  R(t)  can  be  expressed  in  terms  of  the  B 

and  V  autocorrelations  and  the  variances  of  X  and  X'  . 

o  o 

The  maximum  likelihood  estimate  of  x(T)  would  be  obtained  by  maximizing 

p  (U  x)  with  respect  to  x  ,  ,  x  .  and  then  taking  x(T)  =  lim  ix(T) 

r  n°  A  n_i  nvc^i 

where,  for  each  n  .  x(T) I  ■*  x  ,  . 

’  ’  L  ; n  n-1 

This  procedure  appears  to  be  of  only  academic  interest  in  the  general 
case.  Section  V  shows  how  the  problem  of  optimum  position  estimation  can 
be  reduced  to  a  standard  linear  least-squares  filtering  problem.  The  reason 
for  emphasizing  the  maximum  likelihood  approach  at  this  point  is  simply  that 
this  is  the  approach  actually  used  in  obtaining  results  for  Cases  1  and  2„ 
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III.  TREATMENT  OF  CASE  1:  0B(a,t)  “  Y 


A '  m  IVATIOH  0?  OPTIMUM  POSITION  ESTIMATE 

In  this  case,  all  sources  of  error  are  present  except  for  the  non¬ 
bias  component  of  accelerometer  noise. 

Mow,  since  x  (t)  +  fl?  x(t)  “  b(t)  , 

(III.1)  x(t)  «  ~  j  b  (tO  sin/i(t  -D  dr+  x  cos  fit  *  x* 

lit  12 

Since  in  the  present  case,  C>a(t)  *  B(0)  -  b(0), 

O 


(ill *2)  B(t)  -  b(t)  +i^(t)  -  b(t)  +  Bo  -  lxQ 

Let 

J 

rt 

(III. 3)  fit)  -  i  1  B(r)  gin  flit  -  ~)  dr 

j  I'o 


Then,  putting  (III. 2)  into  (III.l), 


(III.U)  x(t)  ■  -J/(t)  +  x  cos  1  It  +  -£  sin  -  it  +  - ^ 

H  0  il  ff 


?',x  -B 

0  0  (1  -  cos  1 it) 


ris  X  —  B 

v(t)  -  x  *  ( t )  “  f'it)  -  xq  fi  sin  fi  t  +  x^  cos  fit  +  — sin  fit 


J*  Z. 


From  (Ill.li)  we  can  say  that 


✓V 


(in. 5)  x(t;T)  »  fit)  *  x(OjT)  cos  fit  *  sin  fit 


A 

<tx(0;T)-3 


IL 


.  -'\ 

(1  -  cos  1  <_t) 
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(III. 5)  -  cont'd) 

v(t;T)  ■  ^(tjT)  =  \f/'{  t)  -  x(0;T)i  )sin  fit  +  P(0;T)  cos  I  it 


4x(0jT)-B 

__2 


Sin 


.Qt 


Applying  the  maximum  likelihood  method  as  explained  in  Section  II, 

one  first  imagines  the  velocity  dial  readings  to  be  given  just  at  times 

tQ  “  0>  ...  ,  t  ^  “  T.  One  must  then  minimize,  with  respect  to'x'(OjT), 

/n  A 

x’(OjT),  and  JCx(OjT)  the  quantity 


(III .6) 


JS.tL_ 

\  <  A 


(V.  -  Q)Ci,  -  i]i  •  ± 


'  »  ,  \  v  l  ■JV«/\VJ  A  I  T  — —  a 

i,j”0  i;3  1  1  i  i  r3  _  o 


X'  -  x(OjT) 


2 

r  /s  l 

♦  i 

Yl4 

X  <  -  x«(0jT) 
0 

—  -J 

+  j_ 
Y1 

V> 

Bq  •*  ir.x(0}T) 

where  (  3  )  a 

1 J 


'‘V^i’  V 


-1 


(i.e„  the  matrix  inverse) 


Qj,  is,  apart  from  a  constant  independent  of  x'(OjT),  x^(0;T)  and 

j  A 

U  I  x(OjT),  x'(OjT),  <Cx(0;T) 


j 

x(OjT),  equal  to  -In  p 


For  convenience  let  x(0jT)  »  z  j  x^(Ojl)  »  z  j  ^x(OjT)  -  B  «=  z  .  On 

•t.  c  0 

differentiating  Qn  with  respect  to  z1,  z2,  z,  (after  substituting  (III. 5)  into 

(ill. 6))  and  setting  the  derivatives  equal  to  zero,  three  equations  are  obtained: 


1  ^7 

(i)  z]_  (An  +  H  -  z2  ^2  "  z3  T  "  D1 

i  L 


(n)  -  z1  A12  +  z2  (A22  +  ♦  z3  -  D2 

U  II 


h2 


(iii)  zl  *  z2  ^12  Z3V“?  *  ’  D1 


+  i.  X 

y3  0 


(III. 7) 


where 


(i) 
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(ii) 


12 


JQrl_ 

■rfi 


-  5  .  .  sin  1  ]t.  cos  -Qt 

i, jHo  x3  1 


'J 


(III. 8) 


n=l_. 


(iii)  A22  -  2____  j±.  cos  .Qti  cog 


1jJ  =  0 


(iv)  \m--rLr-K3<\-f{)* 

i,j=o  -j  i  i  v 


r3  o 


(v) 


^5L_. 

32  "  L~.  -  "  Yi' )  COS  fit.  +  JL  xt 

T7j“o  j  a  '  1  3  ri  0 


The  maximum  likelihood  estimates  x(OjT),  x>(0;T),  ^x(0;T)  -  Bo  when 

the  velocity  dial  readings  are  given  just  at  times  t.  are  obtained  by 
s0lving  Eqs,  (III.7)  for  ^  ^ 

Is  then  obtained  by  putting  the  resulting  values  into  Eq.  (III^J. 

an oi ties  A^,  A^,  A 22 ’  °1»  °2  as  well  as  the  resulting  estim¬ 
ators  depend,  of  course,  on  n.  This  has  been  partially  suppressed  in  the 
notation  for  convenience  in  writing  toe  above  equations.  Eventually  we 

will  evaluate  the  limiting  values  of  these  quantities  as  n~^,  for  a 
particular  j*?(s,t). 

The  limiting  values  of  all  quantities  will  be  denoted  by  the  same 
symbols  as  have  been  used  in  the  above  equations.  It  will  always  be  clear 
from  toe  context  how  to  interpret  the  symbols.  The  limiting  values  of  the 
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A  /\ 

estimates  x(OjT),  x  1  ( 0 j T ) ,  ?>x(OjT)  -  Bo  are  (see  the  remark  in  Section  II) 
the  best  estimates  based  on  the  dial  readings  for  all  t,  0  ®  t- ■”  T.  The 

best  estimate  for  x(t),  for  any  t  in  0  ■  t  "  T,  is  gotten  by  putting  the 

A  /\  <7'  * 

limiting  values  of  xtO;T),  x'(OjT),  and  cAx(0;T)  into  (ill. 5). 


Adding  (iii)  to  (i)  in  (in. 7)  gives 


(III. 9) 


Z~  * 

3 


n'W 

— PX- 

r^L 


or 


a\ 


A, 


£x(0;T)  -  3 


Yi(Xo-x(0jT)) 

o  - 72 - 

r3n 


Substituting  (ill. 9)  into  (HI. 7)  gives 


(i)  i. 


(III.10) 


*11 


1  + 


r3n“/  '3 


AliYl 

“2  ^12  "  D1  +  Xo 


r3rr 

A12Yl 


(ii)  -  A12  ^1  +  +  z2  (A22  +  "  D2  ■  ^ 


• . ■» '  X 

•7?  ° 


Some  special  cases  are: 


(a) 


J 


Yj  n«°(No  independent  initial  estimate  of  xq) 


Iyx  t 


Best  estimate  means  best  in  the  sense  of  the  optimization  criterion  of  • 
Section  I.D  (p.6),  with  x(T)  replaced  by  x(t;T)  and  x(T)  replaced  by  x(t). 
Cf.  the  footnote  on  p.  9. 
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In  this  case 


(111.10a) 


(i)  ^  -  z2  ijj,  -  Dj 


(il)  ■■  zi  ^  +  z2  (A22  *  Y^)  "  d2 


(b)  ,•  “  0  (No  accelerometer  bias) 


(I II. 10b) 


Iyj  f  o 


(i)  zl  (A^  +  Y*  )  -  z2  A12  -  D1 


(li)  -  A12  +  Zg  (Ag2  +  ^)  °  02 


The  form  of  jilt's, t)  has  not  been  utilized  in  deriving  tlie  above 
equations,  which  ere  therefore  valid  for  any  0y(s,t).  The  form  of  $y(s,t) 
affects,  of  course,  the  actual  values  of  the  quantities  A^,  A^g,  A0g,  D^, 
and  Dg.  Appendix  I  shows  how  to  evaluate  the  limits  of  these  quantities  as 
n— >oo,  for  0v(s,t)  of  the  form  r2  +  Pg  e  Is  ^a2  . 

Tlie  limits  of  these  quantities  for  other  0y(s,t)  can  be  evaluated  by 
the  method  of  Appendix  I— and  thus  explicit  solutions  for  the  optimum 
position  estimate  can  be  obtained--provided  one  can  solve  the  integral 
equation  (AI.3)  with  the  appropriate  0v(s,t)  as  kernel. 
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B.  EVALUATION  OF  LIMITING  VALUES  AS  n— )o=for  0^.(s,t)  13  Yg  +  |3^e  ^  ^ 

Appendix  I  shews  how  to  evaluate  the  limiting  values  of  bilinear  forms 


such  as  appear  in  (III. 8).  The  limiting  values  for  the  case 


-a„|s-t| 


VS,t)  "  Y2  +  ?2e  3X6  obtained  by  applying  Eq.  (AI.7)  of  Appendix  I, 

according  to  the  following  key: 


Quantity  evaluated 

Values  of  quantities  appearing  in  Eq.  (AI.7) 

- 2 _ 2 _ Y _  u(t)  y(t) 

A11 

a2 

^2 

I 

Y2 

j  -Osin  Dt 

Osin  Ot 

CsJ 

a2 

h 

Y2 

f - - 

i  Isin  Ot 

coa  Ot 

A22 

- - - I 

°2 

h 

Y2 

n 

C03  j  t't 

i 

cos  Ot 

1 

o 

j 

°2  ! 
- - i 

_ 

Y2 

-  V(t)-y«(t)]  j  Osin  A  it 

d0  -  ~  x  *  ! 

L  2  0  i 

°2  |  h  | 

y2 

-  ■  . i„ 

v ( t )  -  y/i  ( t )  |  cos  Ot 

J. 


The  results  of  this  are: 


.2 


(Ill.lli) 


w-  *  4  V  —  i -A 


“u  \  ~ 


2.  u 


/  i»0  \  a2  / 


sin  2  1  iT 


+  siO  _Qt 


2 

j^-(l-cos  U  It)  +  Sin Ot 


2  +  o„  T  + 

2  Yc 
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(I II *1111) 


*12 "  (f2 1  ro  (J1 "  ^rjsin2  '^T  + 1 oinJ  flT 


•i:-  f 

|y 

§ 


p 

i 


JS 

§ 

!■ 

?K* 

w 

IP 


jj(1  -  cos  n T)  +  sin  Or 

-pp  sin  fll  +  1 

IN 

+  cos  i  It  i 

2po^ 

2  +  a0T  +  — - 

2  r0 

J 

(lll.lliii) 


i  a2T  /  o2  \  °2  /  n? 


I  sin  2  -i  -if  +  1  +  'cos2-Qt 


sin  -Ol  +  1  +  cos  Ql 


2  +  V  + 


2p? 


(Ill.lliv) 


~2 


a2  /  0 


V(t)  -  i  sin  i it 


dt 


v(o)  ♦  (sin  fi T  +  3-i  ccs  XIt) 


7(T)  -Ij/'iT) 


t 


r°: 


in 


%  (1  -  cos  i7t)  +  sin  lIt 


a2  (T  V(TJ  -  l!/( T))  +  V(o)  +  V(T)  -*^(T)  j 


2  +  a,  T  + 


W 
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jgj* 

f? 


ft: 

Sf1? 

I 

I 

I  H 

■'  I- 

ft 

I 

I 

1 


Si 

I 

r 


t 


X 


(Ill.llv) 


D 


<» 

2 


o 


rh 


v(t)  -  y/'(t)  |  cos  _Qt  dt 
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+  V(o) 


jrj  3  In  I'iJ 


tJ  +  1  +  cos 


aii  b  (■ 


TV(T)  -  tMt))*  Ho)  + 


V(T)  -  1y^(T) 


2*“2T*r. 


»h«re  V(T)  ■  ij  J(t)  dt 

'O 

Toe  estimators  x(0jT),  x '(OjT), &(0,T)  -  Bq  approach  limits  which  are 
of  course,  the  solutions  to  Eqs.  (III. 9)  and  (ill. 10)  where  the  various 
quantities  on  the  left  side  now  are  given  by  (III.11), 


C.  .ERROR  STATISTICS: 

faking  mean  values  on  both  sides  of  Eq.  (Ill.io)  yields 

(in. 12) 


+  x ' 
0 
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Applying  these  together  with  Eq.  (III. 9)  gives 


(III. 13) 


°  X  J  S.  “X1  J  2,  °  0 

1  0*2  o  5  3 


Thus  the  estimates  of  xq,  x£,  and  are  unbiased;  also,  therefore, 
the  estimates  of  x(t)  for  any  t,  0  »  t  •  T,  (The  property  of  being  unbiased 
refers,  of  course,  to  ensemble  averages— i.e.,  averages  over  a  large  number 
of  flights.  For  any  given  realization— i.e. ,  single  flight— the  instrument 
biases  may  cause  bias,  in  the  estimate  of  x(t).) 

Now,  for  the  sake  of  brevity,  put 


(i)  Kn”An(y^J^ 


(ii)  K.i2  -  A12 


(iii)  K21UA12  1+ 


r/f; 


(iv)  K22  =  A22  +  ^ 


(III. 110 


1  Y-j 

<T)  b-V-rr 
r3a 

■^*1  o  T*i 

(vi)  E  -  D0  -  -±L-±  x 
2  2  0 


X  - 


^12  Tl 
Do  "  K 

Y 7f 


and  put 


(in. 15) 


/V 


1  h  -  Z1  n  x(a’T>  “  xc 


“  z2  ~  z2  ”  (0;T)  -  x' 
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We  know  that  ^  =  0.  Subtracting  (III, lit)  from  (ill. 12) 

gives 

^  K11  ~  ICL2  "  ^1 

(in. 16) 

(ii)  ~K21  ^  *  K22  ^  ”  E2 


Now,  subtracting  (lll.li)  from  (in. 5)  gives  (for  0  »  t  =  T) 


(III. 17) 


£.A(t;T)  »  JP1  cos  fit  +  J°2  +  -_LL^-5_2  (i  -  cos  fit) 

-»  1  _Q_ 


In  view  of  (ill. 9),  this  may  be  written 


(III .17*)  O^(tjT)  -  J°x 


COS 


Y 

fit  — (i  -  cos  nt) 


'3-' 


7? 


\  u 


+  1°  3in  fit 

w  2  a 


+  c, 


\J:A- 


- — j—  (l  •*  cos  u  it) 


i 

+  — w  (1  -  COS  fit) 

f' ■ 

U  i. 


Thus . 


,  to  obtain  ^/*-(tjT)>  it  suffices  to  know  7°  2 , 7l2 
X  x  2 


7° , 7^  fy  ,  Cg  ,  6g  .  From  (III.  16),  it  is 


seen  that  it  suffices  to  know 


E1  '  E2  »  E1E2  »  E1^C  »  E2°X_  »  V'B  » 


r 

-  6 
“2  B  * 
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?  T~ 


n-1 

5  *  -il  }  5.  . 

1  i7F°  1J 


sin  1  k 
GV(ti}  “  r[ 


sin  fit. 


1 


'  '  *11  Yi 


k*  !  i 


X7 fjh 


I 


and 


nrl_ 

i. 


T77a° 


ij 


-  san  i  ;.:t,  , 

(-W 


COS  -»  it . 


1 6  .  >  ri  f . 

r,  x*  X  'Jx 
i*  O  r  3rf  O 


or 


I  •('} 

i  r 

t  ■ 


(i)  E^  =  -fl} _ _  J.  .  Lv  (t. )  sin  fit. 

X  i,j-o  iJ  V  1  0 


(III. 18) 


a2  B  H  7T/r} 


n-1 


(ii)  E2  “J~~ ?)jLJ  6V  (t.)  cosily 


\2  f  f.2  r'  r 


xtv  .i  a 

o 


n?  B  TcF  "‘o  n,  "si 


The  evaluation  of  the  quantities  E^,  etc.,  can  be  illustrated 
by  the  evaluation  of  E^2: 
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V  -  A  L‘ 


ijkx 


i  ~ij  kj? 


(ti}  CV  (V  sin  fit .  sin  -Qty 


+  in  fi  +  1  f,  +  An  Y1  Y  f  2 

if  ^  r7  v  7f  i  S 


Using  the  fact  that  (t^  (t^)  »  0y  (t^  t.)  and 

f  1  .  i  rr  7T 

L _ '  ik  Jk((  "  811(3  a-^so  *  Y]_j  Oj  *>  one  can 

easily  evaluate  this  expression.  The  result,  'together  with 
the  results  for  the  other  necessary  evaluations,  is: 


(III. 20)  (1)  B?  -  / 1  *  kili  ' 

1  \  IF 


i  ’  ■'Si  yi 
a.-  +  •—  l  i  +  JA..  A  , 

n  t3  \  7F i 


(li)  E!:  “  A„  +  —  +  / 1  +  H 


2  22 


^  jt 


\  YVfy 


In  \  Z'  Tn  K  F 


(ili)  H1E2B-A12  1  +  -~n)ll  + 

V  T-J  )  !  \ 


1  11 

W  V* '  7T  / 


/ 


An  Yi ' 


(iv)  E,  6  -  (  1  +  - r 

1  Xo  V  /?  J 


(v)  E  (/  « 

2  xo  ~^r 


I  c  f  .  ^  Y1 
ivi)  ^  03  ■=  - r 

J  A 
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The  equations  derived  in  this  section  for  mean  values  and  variances 
of  the  various  estimators  hold  in  the  limit  as  n ->  oo  (with ,  of  course,  the 
quantities  ,  etc.,  taking  on  their  limiting  values  as  n-*  °°). 

All  the  relations  (111,12)  -  (III. 20)  hold  for  any  autocorrelation 

function  0y(s,t).  As  noted  above  (p.  17)  the  form  of  0y(s,t)  determines 

the  actual  values  of  A^,  etc, 

jTTT 

The  calculation  of  O^tjT)  is  now  a  straightforward,  though  very  cumber - 
some,  application  oi  the  results  given  by  Eos,  (III. 16),  (III, 17'),  and 
(ill. 20).  3ecause  of  the  cumbersome  nature  of  the  resulting  formula,  it  is 
not  thought  worthwhile  actually  to  write  it  out  here.  However,  certain  con¬ 
clusions  can  be  drawn  fairly  easily.  For  example,  one  may  use  these  equations 
together  with  Eos.  (III. 11)  to  calculate  the  behavior  of  the  estimator  vari¬ 


ances  as  T~>oo,  for  fL(s,t)  “  To  + 


-a2js-t| 


The  results  are:  the  error  variance,  for  sufficiently  large  T,  approaches 
t1Y3 

the  quantity - ; — r  .  (This  checks  a  more  general  result  obtained  in 

/->U 

Y1  +  y3‘  ± 

Section  VI,  namely  that  this  quantity  is  the  asymptotic  error  variance  in 

( 

Case  1  for  .any  such  that  C^(t)  is  stationary  and  has  a  power  spectrum 
containing  no  delta  function  at  frequency  d  l ,)  For  y-j_  or  Yj  ”  0,  the  error 

x  -s!2-tl 

variance  goes  to  zero  as  ^  for  (fly(s,t)  c  +  p.e 


D,  SUMMARY  OF  RESULTS,  CASE  1 

Tne  best  estimate  of  position  at  time  T,  based  on  all  information  received 
up  to  time  T,  is 


r>- 


c6  x  ( 0 ;  T )  -  B 


(III. 21)  x(T)  -  B(V)  -  -  +  - — * - -  (1-cos.,  t) 

!  ,  —  ^ 


/V 


/\ 


in 


x(0;T)  cos  -  AT  +  x'(D;7)  — 


CONFIDENTIAL 


0 


lunriutNiiAL  ™ 

3-25-pii 

a,  /\  44  ”26- 

where  x(OjT),  z'(OjT),  Sq  -  ^x(OjT)  are  respectively  the  best  estimates 

Ox  initial  position,  initial  velocity,  and  accelerometer  bias,  based  on 

all  information  up  to  time  T„  x(OjT),  x'(0;T),  and  <£x(0;T)  -  Bo  are 

solutions  of  a  set  of  simultaneous  linear  equations  ' Eqs .  (m,9)  (m  10) 

A  /\  ''N  ' 

with  z  m  x(0jT)j  x  •  x'(OjT);  and  z,  -  Xx(OiT)  -  B  ). 

This  holds  for  any  velocity  noise  autocorrelation  function  0  (8, t). 

The  function  0v(s,t)  affects  the  values  of  the  coefficients  of  the  simul¬ 
taneous  equations.  For  0^(S,t)  -  Tg  ♦  ^e“d2  M  ,  ^  ^  ^ 

found  from  Eqs.  (ill. 11). 

In  other  words,  x(T)  is  found  by  solving  the  accelerometer  equation 
(Eq.  (1.1))  with  B(t),  the  accelerometer  dial  reading,  as  driving  function; 
subtracting  a  term  which  represents  the  best  estimate  of  the  effect  of 
accelerometer  bias;  and  using  as  initial  conditions  the  best  estimates  (up 
to  time  T)  of  initial  position  and  velocity. 

As  ^lapsed  time  Throws  larger,  the  computer  continuously  re-evaluates 
x(0;T),  x'(OjT),  and  Xx(0;T)  -  Bo.  (For  0  •  t<T,  one  can  get  the  best 
estimate  of  what  the  position  was  at  time  t  from  Eq.  (ni.5).  Thus,  for 

any  fixed  t<T,  the  computer  can  be  used  to  get  a  continually  better  esti¬ 
mate  of  x(t)  a3  T  grows  larger.) 

The  estimates  x(0;T),  O(0;T),  and  £(0;T)  -  3q  have  eixors  which  have 
sero  mean  (over  a  large  number  of  flights),  are  linear  functions  of  the  dial 
errors,  and  do  not  depend  on  the  path.  Thus 

(ill .22)  2(T)  -  x(T)  ♦  ^(T) 

where  t^(T)  has  zero  mean,  is  a  linear  function  of  the  dial  errors,  and  is 
independent  of  the  path.  The  variance  of  f^(T)  can  be  found  from  Eqs.  (ill. 16), 
(III. 17'),  and  (ill. 20) 
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c 


F 


L 


IV.  TREATMENT  OF  CASE  2:  “  Y2 

A .  DERI VATIOM  OF  OPTIMUM  POSITION  ESTIMATE 

In  this  case,  all  sources  of  error  are  present  except  for  the  non¬ 
bias  coinponent  of  accelerometer  noise. 

Since  x'(t)  "  v(t),  we  have 

(IV.1)  x(t)  °  )  ▼  CC)  dt  +  XQ 

'o 

Also  in  this  case  6y(t)  a  V(0)  -  v(0),  so  that 

(IV.2)  V(t)  -  v(t)  +  Vq  - 

Let 

(IV. 3)  W(t)  “  j  V(C)  dK 

v 

then,  putting  (IV, 2)  into  (IV. 1), 

(IV. U)  *U)  ■  W(t)  +  xQ  +  (xj  -  7q)  t 

and  ^x(t)  “  aSW(t)  +  -0.  xq  +  (xj  -  Vq)  fl  t 
Thei'efore 

(IV. 5)  x(t;T)  »  W(t)  +  x(OjT)  +  (x^(0;T)  -  Vq)  t 

A  2  /\  „  2 

£x(t;T)  -  cTW(t)  +  -ft  x(0;T)  +  (x  1  ( 0 ;T)  -  VQ)  12  t 
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*11 


-  n  i 


n-L 


T'ij 


(IV. 8) 


*12 


n£L_ 


1jJ"0 


/  A 


22 


ft£L 


i  J"o 


T'«  ‘i  kj 


2 

Di  ■  A- 1~-  \i  (Bi  ■  '  V  * 


i,J=o 


--  \  *  X'-V 

Dp  -  U  /  ru  .  (B  -  cCw  )  t.  ♦  -2—2 

2  t^=o  1  1  J  *0, 


\ 
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The  best  estimate  of  x(t),  0  =  t  “  T,  when  the  accelerometer  dial 
readings  are  given  just  at  times  is  now  obtained  by  putting  the  solu¬ 
tions  to  (IV, 7)  into  (IV. 5).  As  in  Case  1,  we  will  evaluate  the  limiting 
values  of  the  solutions  i  ,  z of  (IV.7) ;  the  best  estimate  of  x(t),  0  «  t  £  T, 
will  then  be  obtained  by  putting  these  limiting  values  into  (IV.5). 

The  above  equations  are  valid  for  any  0  (s,t),  The  method  of  Appendix  I 
can  be  applied  tc  evaluate  the  limits  of  A^,  A]2,  k22,  D1,  and  D2  for 

-Is-t la. 

0g(s»t)  *  Yj  +  Pxe  1  1  i. 

The  limits  of  these  quantities  for  other  ^ ( s , t )  can  be  evaluated  by 
the  method  of  Appendix  I  provided  one  can  solve  the  integral  Eq,  (AI.3) 
with  the  appropriate  0_(s,t)  as  kernel. 
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+>  leg 


(lV.9iv) 
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1  “  2^  “lT  +  Bo 


?>W  +  B(T) 


^W(T)  1 


(!V.9v) 


2  f  i'T 

°2  ’  °1  jo  1  B(t)  -  dt  -  ~  (Bq  -  ^WQ)  +  (T  +  i.)  B(T)  - 


f  dM- 

-4 

r  . 


£  a^T  (B  -  c-SW)  +  Bq  -  cAW  +  B(T) 


X  -V 
o  o 


where 


B  -  c^W  n  -i  J  ["  B(t)  -  ebW(t)  dt 


C.  ERROR  STATISTICS: 

Taking  mean  values  on  both  sides  of  (IV. 7)  gives 


(i)  %  <An  *  7j  *  z2  *12 '  ‘'li  *  ^  xo 


(IV.10) 


(ii)  zx  A12  +  z2  (A22  +  -L  +  -£.)  =  ^2  xc 


Applying  these  equations  gives 


Z1  “  Xo  5  Z2 
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(IV.12) 


1  “  *1  -  **  x(0;T)  - 


A 


2  =  z2  "  "2  "  x 1 (°JT)  -  V 


(IV. 13) 


gives 


E2  "  °2  -  °2 


(IV.lil) 


(i)  yi  <*U  *  *  “'s  *12  ■  E1 

(u)  j\  Au .  3;  (a22  ♦  i  .  i)  ■  Ej 


Now,  subtracting  (iV.ij)  from  (IV.?)  gives 


(17.15) 


■^(t;T)  ■  (x(0;T)  -  x)  +  (x^(OjT)  -  x»)  t 


^1  +  (v/2  +  V  1 


?  £ 


Therefore,  to  obtain  CaUjT),  one  needs  to  knew 


^ \  »  w^2  }  '^1  ^2  '  5  '"^2  W  *  Fro::i  Ecj*  (1^*^)  it  iu 

evident  that  it  suffices  to  know  ,  Eg2,  ^  !y  ,  Eg  8y  . 
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Now 


nrl- 


^  EfriZ _ n,. 

i>j-o  2 
n-1 


-B(ti}  “  ^Zti  ~v 


'-X 


(IV.16) 


S  *!2  *  TT 


2  ^  G 

(ii)  E2  “  ^  Z _  1\1  6B(tt)  t  -  (A  +  ~)  £  +  _J> 

i,j-o  8  1  J  22  V  v  Ti, 


Evaluating  the  necessary  quantities  in  a  manner  similar  to  that  used, 
in  Case  1,  one  obtains: 


U)  ^  i  *  r2  ^ 


(il)  522  "  *22  *  ^  *  V  S?  *  k  ' 


(iv. 17)  (iii)  »  A12 


Yli  '2^22  ru 


1  +  — -  +  y  £ 

Tjj  T2  *22 


l 


(iv)  E,  Otr  *■  - 


1  W  "  ■  r2  A12 


W  s2^“-r2(A22.i 

V  V 


All  the  relations  (IV.10)-(IV.17)  hold  for  any  autocorrelation  function 
The  form  of  ^(s,m)  determines  the  actual  values  of  etc. 

All  these  relations  hold  in  the  limit  as  n->oc. 
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The  error  variance  can  be  caJ.culated  in  a  straightforward  mannei  using 
the  above  equations.  As  in  Case  1,  the  resulting  formula  will  not  be  stated 

because  of  its  cumbersome  nature. 

Applying  the  above  equations  together  with  (TV. 9)  for 
-a, |s-t| 

0  (s,t)  “  y  +  {Le  ^  j  one  obtains  the  result  that  the  asymptotic  error 
B  9  13. 

Y  y 

variance  for  large  T  is  - — r  •  (ibis  checks  a  more  general  i  esult 

obtained  in  Section  VI,  namely,  that  this  quantity  is  the  asymptotic  error 

variance  in  Case  2  for  any  such  that  <Sg(t)  is  stationary. ) 

For  Yi  or  y-j  d  Oj  the  error  variance  goes  to  zero  as  ^  f°r 
1  J  -cujs-tl 

i^(s,t)  «  yx  *  Pxe  '• 

D0  SUMMARY  OF  RESULTS ,  CASE  2 

The  best  estimate  of  position  at  time  T,  based  on  all  information  received 
up  to  time  T,  is 

(IV.18)  x(T)  -  |  V(T)dt*  x(0;T)  +  (£(OjT)  -  Vo)T 

Jo 

where  x(0;T)  and  Vq  -  jf'(0$T)  are  respectively  the  best  estimates  of  initial 

position  and  of  velocity  dial  bias,  based  on  all  information  up  to  time  T» 

x(0$T)  and  x»(0;T)  -  Vq  are  solutions  of  a  pair  of  linear  simultaneous  equations  (Eqs. 

(IV. 7),  with  z1  *  x(0;T);  “  x^(0;T)  -  VQ). 

This  holds  for  any  accelerometer  noise  autocorrelation  function  0g(s,t). 

Hie  function  fi(D(s,t)  affects  the  values  of  the  coefficients  of  the  simultaneous 
B  -ojs-t) 

equations.  For  0  (s,t)  "  Yx  +  Pxe  ,  these  values  can  be  found  from 

Eq.  (IV.9). 

Thus,  x(T)  is  formed  by  integrating  the  velocity  dial  reading  V(t); 
subtracting  a  term  which  represents  the  best  estimate  of  the  effect  of 
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velocity  dial  bias}  and  using  as  initial  position  the  best  estimate  of 

x  based  on  all  information  received  up  to  time  T. 
o 

The  estimates  x(0|T)  and  x^(0;T)  have  errors  which  have  zero  mean,  are 
linear  functions  of  the  instrument  errors,  and  do  not  depend  on  the  true  path. 
The  same  things  therefore  are  true  of  the  error 


(17.19) 


x(T)  ■  x(T)  +  8a(t) 


where  ua^  has  mean  zero,  is 
is  independent  of  the  true  path 
from  Eqs.  ( IV .lU ) ,  (IV. IS),  and 


a  linear  function  of  instrument  errors,  and 
x(t).  The  variance  of  ^(T)  can  be  found 
(IV. 17). 
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I  \  “  R(0)  ♦  s(o) 

I 

j 

!  »  R'(o) 

(V.2)  J 

\ 

V(t)  -  R'(t)  +  S'(t) 

I  ^2 

\  B(t)  -  R"(t)  ♦  -L  R(t) 

Thus  the  problem  of  obtaining  an  estimate  x(t)  for  x(t)  from  a 
knowledge  of  B(t),  V(t),  XqJ  X^  is  completely  equivalent  to  the  problem 
of  obtaining  x  t)  from  a  knowledge  of  any  two  of  the  three  dials  h(z) s 
C(t),  S(t). 
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Note  that  Q(t)  and  R(t)  each  are  estimates  of  x(t);  S(t)  is  the 
difference  between  these  two  estimates  of  x(t). 

The  advantage  of  the  reduction  from  four  dials  to  two  is  that  it  is 
simpler  to  visualise  intuitively  what  the  optimum  process  amounts  to  if 
one  considers  the  input  dials  to  be,  say,  R(t)  and  S(t),  rather  than  the 
original  four  dials.  Also,  for  the  purpose  of  calculating  asymptotic  errors 
(see  Section  VI),  the  reduction  turns  out  to  be  mathematically  convenient. 

If  we  define 


4q(t)  "  U  +  1  C  (r)  d'C 

o  B 


(V.3)  ^£R(t)  B  cos 


/•  , 

°X'  - 

fit  +  sin  fit  ;  6a(r)  sin  flit  -r)  <n 

o 


Lit) »  c  (t)  -  £_(t) 


vts  are  led  to  the  result 


(V.h) 


S' 

/ 

!  R(t)  °  x(t)  +  L(t) 

ft 


(  Q(t)  »  x(t)  +  £Q(t) 


|s(t)  -  £s(t)  -  6Q(t)  -  iRit) 

V 


The  problem  assumes  a  particularly  convenient  form  if  one  considers 


the  pair  of  oials  [f?(t),  S(t)j  or  the  pair  j  Q(t),  S(t)j.  For  expository 
purposes ,  we  will  consider  the  pair  f.(t),  S( t)  .  These  will  be  referred 
to  as  eigen dials. 


1 
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From  (y.lt)  it  is  apparent  that  an  estimate  of  x(T)  can  be  obtained 
by  estimating  £r(T)  on  the  basis  of  the  readings  R(t),  S(t)  for  0  »  fc  -  'f , 
and  then  subtracting  this  from  R(T).  In  other  words,  one  can  define  an 
estimate  of  x(T)  by 

A 

(V.5)  x( T)  -  R(T)  -  6r(T) 

A 

where  £r(T)  denotes  the  estimate  of  6  (t)  based  on  R(t),  S(t)  for  0  -  t  -  T. 
From  this  it  follows  that  (subtracting  (V.h)  from  (V.5)) 

A 

(^*6)  Error  in  x(T)  B  Error  in  6  (T) 

R 


It  is  also  true  that  all  possible  estimates  of  x(T)  in  terms  of  R(t), 

S(t.)  for  0  -  t  -  T  are  of  the  form  (V.5)  for  suitably  defined  Gd(T),  since, 

- 

given  an  estimate  x(T) ,  one  can  define  the  estimate  of  <bR(T)  by  £r(T)  =  R(T)  -  x(T). 

We  can  conclude  that  any  process  of  estimation  of  x(T)  on  the  basis  of  the 
dial  readings  for  0  -  t  =  T  is  equivalent  to  first  estimating  £_(T)  on  the  basis 
of  the  eigendial  readings  for  0  °  t  =  T  and  then  obtaining  x(T)  by  (V.5).* 


B.  REDUCTION  TO  A  FILTERING  PROBLEM 

According  to  what  was  said  in  the  introduction,  one  desires  to  obtain  a 
linear  estimate  of  x(T)  such  that  ^x(T)  -  x(T)J  is  a  miniumim,  and  this  esti¬ 
mate  must  not  utilize  any  a-priori  information  about  the  possible  paths  x(t). 
In  view  of  (V»5)  and  (V.6)  the  first  requirement  means  that  one  desires  to 
lind  an  estimate  *-r(T)  of  £r(T),  formed  by  a  linear  operation  on  R(t),  S(t) 
for  t  T,  for  which  j^„(T)  -  CR(T)J2  is  a  minimum.  The  second  require¬ 
ment  obviously  implies  that  R(t)  yields  no  information  about  £R(t)  -  i.e. 


This  is  closely  related  to  a  more  specialized  result  of  J.  H.  Lanina- 
see  Ref.  ?. 
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that  the  estimate  of  C  (T)  joust  be  based  on  S(t)  alone.  In  short,  the  problem 

K 

amounts  to  one  of  linear  least  squares  filtering  of 
6R(t)  from  S(t)  “  ^q(t)  -  CR(t).* 

C.  ALTERHATO  DSFIMiTIOH  OF  EIGSHDIAI5 

For  some  purposes  it  may  be  more  convenient  to  adapt  Q(t),  S(t)  as 
eigendials.  In  this  case,  one  must  find  a  linear  least-squares  estimate 
of  60(T)  base  on  S(t),  0  ■  t  “  T.  The  equation  for  the  position  estimate 
is  now 


(V.7) 


/\ 

x(T)  -  Q(T)  -  6q(T) 


if 

1 


(V.9) 


Also 


A 

Furor  in  x(T)  b  Error  in  Cq(T) 


For  other  purposes  it  might  be  best  to  use  R(t)  and  Q(t)  as  eigendials. 
This  might  be  the  case  if  one  wished  to  apply  the  maximum  likelihood  method, 
(See  p.  12) 

It  should  be  understood,  of  course,  that  these  various  possibilities 
are  merely  tantamount  to  alternative  mathematical  notations,  and  have  no 
effect  on  the  final  answer. 


% 


If  it  were  desired  to  take  advantage  of  a-priori  information  about 
x(t)  then  P (I )  would  have  to  be  based  on  both  R(t)  arid  S(t). 
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VI.  ASYMPTOTIC  VARIANCF,  OF  THE  POSITION  ESTIMATE 

A.  BACKGROUND 

For  some  guidance  applications  it  is  of  interest  to  consider  the  variance 
of  x(T)  for  values  of  T  large  compared  to  l/jQ  .  This  section  is  devoted  to  a 
discussion  of  this  subject. 

It  is  convenient  to  distinguish  between  the  bias  and  non-bias  components 
of  the  dial  errors: 

Let 


(VI .1) 


C 

Ky  B  Velocity  dial  bias;  Uy(t)  *  Lly(t)  - 


Kg  “  Accelerometer  bias;  £fl(t)  “  C0(t)  -  Kg 


According  to  the  preceding  section,  the  estimation  of  x(T)  is  equiva¬ 
lent  to  estimating  L_(T)  by  operating  on  S(t),  0  =  t  ®  T,  and  then  setting 

A.  W  A 

x(T)  »  Q(T)  -  Cq(T).  The  error  in  x(T)  then  equals  the  error  in  &q(T). 
The  asymptotic  variance  of  x(T)  will  depend  on  the  accuracy  with  which 

f  t 


(VI.  2) 


r\j 

Gn(t)  B  +  Kyt  +  ;  ^.(tO  df 


can  be  filtered-  from 


(vi. 3)  s(t)  -|£x  -  •—  j (l-cos  fit)  *  Kyt  6vw  dr 


L 


X'  -  rt 


a 


-  sin  .fit  -  —  f  ^D(tr)  sin  ri(t-t)  dt 


n)0  B 
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It  will  be  assumed  throughout  this  section  that  &y(t)  and  ^(t)  are 
stationary  processes,  possessing  autocorrelation  functions a 


3.  C/iSE  WHERE  8n(t)  -  0 
B 


A  process  of  operating  on  S(t)  will  be  described,  by  which  all  terms 

of  (VI, 3)  except  for  can  be  extracted  exactly  from  S(t)  as  t-Jo^j 

‘  o  ( 

and  by  which  the  problem  of  estimating  C.-q(T)  can  be  reduced  to  the  trivial 

C 

problem  of  estimating  from  a  knowledge  of  Cy  ~  n  • 

o  o  ff 


(n.k) 


S 1  (t)  “  Ky  +  Cy(t)  +  “  -^jf-O-ain -Qt 


L  .  cos  .fit 
o 


From  (VI.3),  Ky  can  be  obtained  by 


(VI. 5) 


Ky  B  lira  ipi 

v  OO  1 


Also,  assuming  the  spectrum  of  ^y(t)  has  no  delta,  function  at  frequency  -  l9 
one  has,  by  (VI. U), 


(VI. 6)  Ov  ~  “5  "  2  lim  tjt  i  S*(t)  - 

xo  a  *>*>  'J  t 


>in  l  it 


(’/1. 7) 


,  ,T  r 

- 5  .  .  2  lim  i  !  S ' ( t)  -  K  cos  Tat  dt 

a  t-kct;  -  VJ 

0 
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So 

A 

( 

■ 

Oq(1) J 

ffe 

1  ’ 

Sl: 

(VI .11) 
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lim 

T~*oa 


x(l')  -  x(T) 


Y1  r3 


Yi  -  r3  ff 


C.  CASE  WHERE  t  (t)  -  0 


In  this  case 


(VI. 12)  S"(t)  +  fls(t)  -  fl  (i  -\  )  + 


Also 


(VI .13) 


^)  +  Kvru- 

i  ij 


^(t) 


K  -  lia  T 

T-405  1 


S(T) 


and 


t  ?r 


(vi.lli) 


^  -  lim  if 

1  'r. 


•  T 


Xo  Cl  T_>  550 


S"(T) 

a* 


s(t)  -  yt) 


dt 


The  remainder  of  the  discussion  is  the  sane  as  in  the  previous  case. 
Eq.  (VI. 11)  holds  in  this  case  also. 


D.  GENERAL  CASE 

A  ~ 

In  the  general  case  it  is  not  possible  to  separate  f  (t)  »  t  (t)d~ 


exactly  from 


fE(t)  - /  t3m  smCf'-v 


dV  even  as  t — >  o o. 


'o  fl 

Consequently  the  variance  of  x(T)  as  T->  oo  Hill  be  larger  than  in  the  two 
special  cases  discussed  above. 
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Let 


^y(®)  “  spectrum  of  £  (t) 

-r  ^ 

$g(«)  °  spectrum  of  C  (t) 


Since  fy(t)  and  fg(t)  are  not  stationary,  they  do  not  possess  spectra 
in  the  strict  sense.  However,  heuristically  one  should  think  of  f  (t), 
fB(t)  as  having  spectra 


It  seems  reasonable  to  conclude  that  in  the  general  case 


(VI .15) 


lim 

TS'oo 


x(T)  -  x(T) 


n 


Yi  r3 


+  Additional  term 


where  "Additional  term"  depends  on  Jy( »)  and  Jb(cd)  in  such  a  way  that  it 
is  small  if  the  graphs  of  ^(co)  and  ^(<o)  do  not  overlap  much: 


Tig.  2 
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EVALUATION'  OF  LIMITS  OF  BILINEAR  FORMS 


The  problem  at  hand  is  to  evaluate  P  n  lim  P  ,  where 

„  n 

n-400 


RM-1220 

3- 25-Sli 

-1,5- 


(AI.l) 


P”  ‘  b=o  51 


Here  u5  -  u^),  y  -  y(tj)j  t±,  i  -  0,  1#  ...  ,  n  ■  1  are  n  equally- 
spaced  points  in  the  interval  (0,T)  with  t&  =  0,  t  **  T.  u(t)  is  a 
continuous  function  for  0  =  t  »  T,  y(t)  is  an  entire  function  of  t. 


*^o  [\J-  [(Wvtj)]'1 


■at,  -t 

where  0(t. ,t  )  «■  y  +  pe  1  J 

0-  j 


As  pointed  out  in  lief.  1},  pp.  295-296, 


(AI.2) 


f  T>0 

P  “  J  u(t)  w(t)  dt 

-0 


where  w(t)  is  the  solution  of 


(AI.3) 


— « | s-t I 

T  ♦  pe  w(s)  ds 


To  find  the  solution  of  (AI.3),  rewrite  (AI.3)  as 


(AI.lj) 


-a  s-t 

y(t)  -  I  ■  {3e  w(s)  ds 

-0 


where  I  «=  y  J  w(s)  ds 

CONFIDENTIAL 


ss  s  s  e  j 


SjM-1220 

3“25-a 

-4.6“ 


If  we  forget  for  the  moment  that  I  actually  depends  on  the  unknown 
function  w  then  (Al.b)  is  the  same  equation  considered  in  Ref.  b.  The 
solution  is  given  in  Ref.  li*  Eq.  50s 


(AI.5) 


w(t)  ** 


y(t)  -  I  -  4  y"(t) 


J 


2P 


y( 0)  -  1  -  -  y'(o) 


6  (t-o) 


+  1 


y(T)  -  I  +  ~  y'(T) 


6  (t-T) 


Now  integrate  both  sides  of  (AI.5)  over  (-Q,T+0)  and  solve  for  I.  The 
result  is 


(AI.6) 


I  « 


rT 

c j  y(t)  dt  +  y(0)  +  y(T) 

2  '+  ar  +  H 

r 


To  obtain  P,  substitute  (AI.6)  into  (AI.5),  and  then  substitute  the 
resulting  expression  for  w(t)  into  (AI.2).  The  result  is 


n 


(AI.7)  P  =  ^  j  u(t) 


y(t)  -  yn(t) 


dt  +  ^  u(0)  7(0)  -  |  y'(0) 


2? 


u(T) 


y(T)  -  iy'(T) 


1  [a  T  u(T)  +  u(0)  +  u(T)jja  T  y(T)  +  y(0)  +  y(T) 


2  +  c  T  + 

r 


where 


T 


L 


q(T)  “  i  u(t)  dt,  y(T)  *=  i  '  y(t)  dt 
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CASE  3:  0g(s,t)  -  p^e  5  0y(s,t)  -  P2e 

The  method  outlined  in  Section  II  for  dealing  with  this  case  seems 
exceedingly  difficult  to  carry  through.  Therefore,  a  method  more  amenable 
to  analysis  is  used  in  the  following}  this  method  is  still  based  on  the 
maximum  likelihood  principle. 

Consideration  will  be  restricted  to  dial  readings  B(t),  V(t)  which  are 
continuous  functions.  With  the  assumed  autocorrelation  functions,  it  is 
possible  to  do  this. 

The  particular  method  used  below  of  deriving  the  optimum  position 
estimation  process  also  reqaires  that  B(t)  and  V(t)  have  a  certain  number 
of  derivatives.  Now,  with  the  assumed  autocorrelation  functions,  B(t)  and 
v(t)  will  be  nondifferentiable  with  probability  one.  Hence  the  particular 
derivation  given  below  will,  strictly  speaking,  be  valid  for  a  set  of  functions 
B(t),  V(t)  having  probability  zero.  However,  the  final  answer — i.e.  the  method 
of  obtaining  x(T)  from  B(t),  V(t) — involves  operations  on  3(t)  and  V(t)  which 
are  valid  for  any  continuous  3(t),  V(t),  Also,  ary  continuous  function  can  be 
approximated  in  (0,T)  as  closely  as  desired  by  functions  having  any  desired 
number  of  derivatives.  These  facts  are  assumed  to  imply  that  the  final  answer 
is  the  optimum — i.e.  linear  least  squares— estimate  of  position  for  any  contin¬ 
uous  dial  readings  B(t),  V ( t ) .  Tide  amounts  to  assuming  that  if  3j,(t),  V^lt) 
yield  linear  least  squares  estimates  x^(T)  j  if  B(t),  V(t)  yield  a  linear 
least  squares  estimate  x(T);  and  ii  B^Ct),  V^(t)  approach  B(t),  V(t)  in  (0,T) 
in  some  appropriate  sense,  then  x(T)  is  the  limit  of  x^T), 

The  gist  of  the  above  remarks  is  that  various  limiting  processes  will 
be  interchanged  when  such  seems  intuitively  justifiable. 
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As  in  Sections  III  and  IV,  the  maximum  likelihood  method  will  be  applied 

with  Gaussian  statistics.  The  result  its  then  the  linear  least  squares  estimate 

for  any  statistics  having  tine  same  means  and  correlation  functions. 

Consider  the  random  vector  1?  ■  (B  ,  B  , ,  V  ,  ....  V  , .  X  ,  X'). 

o’  '  n-1'  o '  '  n-1'  o’  o  ' 

where  10  B(t,  );  B  V(t^).  The  logarithm  of  the  joint  probability  density 


of  U  is,  apart  from  an  additive  constant,  equal  to  - 


x(t) 


,  where 


(AII.l) 


>\m  L- 
J  i» 

n-1 

V 


«  (vi  -  -  *j>  *  r3  (xo  ■  V2 


*4-^  "u  <Bi  -  J'*t>  <»i  -  S>  *  ^  -  xi)l 


Here 

5  ’ 

ij 

a 

0,  <vy 

L  j 

d  j 

OI 

\  <vy] 

-1 


1-1 


The  limit  oi  this  quantity  as  n-)  c<c  can  be  evaluated  by  the  method 
of  Appendix  I.  The  result  is  (after  suitable  integrations  by  parts) 


(All. 2)  Q 


x(t) 


“l 

f 

/  2  ,  1  ,2 

..  i 

0 

;  U  +  -j  U<  j 

V  °i  J 

dt  +  jp- 

/  T 

zq  J, 

o  1 

(2  +  4  y ) 

\  4  / 

dt'^ 

u2(0)  +  u2(T) 


y2(0)  ♦  y2(T) 


+  —  (X  -  X  )2  +  —  (X  1  -  x')2 
Y  0  o  Y  '  o  o 

3  'h 


whei'e  u(t)  **  B( t)  -  «Ax(t) 
y(t)  -  V(t)  -  x'(t) 
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To  find  the  linear  least  squares  estimate  of  x(t),  0  =»“  t  ■  T,  based  on 
the  dial  readings  up  to  time  T,  one  must  minimize  the  quantity  Q  jc(tjT) 
with  respect  to  the  function  x(t;T). 

This  can  be  done  by  the  usual  Euler  equation  with  free  end  conditions. 
The  result  is: 

The  Euler  equation  satisfied  by  x(t}T)  is,  in  symbolic  notation,* 


(AH. 3) 


Vi 


■j-L-  $(&  -  4)  jB(t)  -  ]  .  A  «t) 


In  addition,  there  are  six  end  conditions  involving  the  values  of 
x(tjT),  B(t),  V(t),  and  their  derivatives  at  t  ®  0  and  t  ®  T. 

Tlie  simplest  way  to  utilize  the  end  conditions  seems  to  be  as  follows: 

The  general  solution  to  (All. 3)  consists  of  certain  integral  transforms 
of  B  and  V,  plus  an  expression  involving  the  initial  values  of  £(t;T)  and 
its  first  five  derivatives,  together  with  solutions  of  the  homogeneous 
equation  corresponding  to  (All, 3).  The  Euler  end  conditions  can  be  used 
to  boil  (his  expression  down  to  one  involving  only  x(OjT),  x'(OjT),  and 
^x(Oji)  -  Bq.  The  function  x(tjT)  thus  obtained  can  be  re-substituted  into 
(All, 2)  and  the  resulting  Q  minimised  with  respect  to  x(0;T),  x>(0;T)  and 
<Kx{ OjT)  -  B  . 


The  final  result  can  be  expressed  as  follows: 


£  ^ 

All  derivatives  of  x(t;T)  referred  to  mean  derivatives  with  respect  to  t 
T  being  regarded  as  fixed.  ’ 


i 
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(AII.li)  *(tjT) 


'^(t  -r)  B(r)  -  K2(t  -r)  v(r)J  dr 


+  x(0jT)  fx(t)  +  x'(OjT)  f2(t) 


<£x(0;T)  -  B  1  f3(t)  +  Xq  fu(t) 


*  K  +  v°  f6(t) 

where  f  (t),  ....  ,  ffi(t)  are  solutions  of  the  homogeneous  equation  corres¬ 
ponding  to  (AII.3).  The  kernels  ^(t),  Kg(t)  can  best  be  given  by  giving 
their  Laplace  transforms i 


(All. 5) 

Laplace  transform  of  K^(t) 


1  ■'  2  + 

*ih 


'  2  2' 


Laplace  transform  of  ^(t) 


1  /  2  2 

=APvP  -*» 


"ib. 


2  ,  n2'' 
p  +  A  t 


2 


2  2 
P  "  <*]_ 


1  2/2  2  \ 

-  ?  (p  "  a2/ 


A 


The  quantities  x(OjT),  x'(OjT)  and  c^x(0;T)  -  Bq  are 
the  solution  to  a  sot  of  linear  simultaneous  equations: 


(All. 6) 


k**l 


A z.  =  D . 

jk  k  j 


j  °  1)2»3 
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wh0re  z±  ~x(0;T)  j  z2  -  £*(0;T)  ;  -  £j(0;T)  -  Bq 

The  quantities  depend  on  T  and  on  the  functions  ...  }  f , 

{as  well  as  a2,  P2)j  Dj  depend  on  these  and  on  the  dial  readings 
XQ,  7 l,  Bit),  V(t),  0  »  t  i  T. 

All  operations  involved  in  finding  x(t;T)  are  valid  for  any  continuous 
functions  B(t)  and  V(t). 
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